Domain wall mobility in nanowires: transverse versus vortex walls 



o 
o 

(N 

> 

o 
O 

m 



> 

(N 
O 

m 
O 



c 

o 
o 



13 



R. Wieser, U. Nowak and K. D. Usadel 
Theoretische Tieftemperaturphysik, Gerhard-Mercator-Universitat Duisburg, 47048 Duisburg, Germany 

(Dated: February 2, 2008) 

The motion of domain walls in ferromagnetic, cylindrical nanowires is investigated numerically by 
solving the Landau-Lifshitz-Gilbert equation for a classical spin model in which energy contributions 
from exchange, crystalline anisotropy, dipole-dipole interaction, and a driving magnetic field are 
considered. Depending on the diameter, either transverse domain walls or vortex walls are found. 
The transverse domain wall is observed for diameters smaller than the exchange length of the given 
material. Here, the system behaves effectively one-dimensional and the domain wall mobility agrees 
with a result derived for a one-dimensional wall by Slonczewski. For low damping the domain wall 
mobility decreases with decreasing damping constant. With increasing diameter, a crossover to 
a vortex wall sets in which enhances the domain wall mobility drastically. For a vortex wall the 
domain wall mobility is described by the Walker-formula, with a domain wall width depending on 
the diameter of the wire. The main difference is the dependence on damping: for a vortex wall the 
domain wall mobility can be drastically increased for small values of the damping constant up to a 
factor of I /a 2 . 

PACS numbers: 75.10.Hk, 75.40.Mg, 75.60.Ch 



Arrays of magnetic nanowires are possible candidates 
for patterned magnetic storage media Q, Q . For these 
nanowires and also for other future magneto-electronic 
devices the understanding of domain wall motion and 
mobility is important for the controlled switching of the 
nanostructure. In a recent experiment the velocity of a 
domain wall in a NiFc/Cu/NiFc trilayer was investigated 
using the GMR effect Q. The measured velocities were 
compared with the Landau-Lifshitz formula for domain 
wall motion 4] . This comparison was used to determine 
the damping constant of the trilayer, a quantity which 
is usually not known a priori. However, several formulas 
for the velocity of a domain wall can be found in the lit- 
erature 0, 0, IE S @ which are derived in different limits 
and all in (quasi) one-dimensional models neglecting the 
possible influence of non-uniform spin structures within 
the domain wall. Thus the question arises in how far 
these formulas are applicable to real three dimensional 
domain structures. To shed some light onto this problem 
we numerically investigate the domain wall mobility in 
nanowires starting from a three dimensional local spin 
model. 

In the following we consider a classical spin model 
with energy contributions from exchange, crystalline 
anisotropy, dipole-dipole interaction, and a driving mag- 
netic field. Such a spin model for the description of mag- 
netic nanostructures 9] can be justified following differ- 
ent lines: on the one hand it is the classical limit of a 
quantum mechanical, localized spin model, on the other 
hand it might be interpreted as the discretized version 
of a micromagnetic continuum model, where the charge 
distribution for a single cell of the discretized lattice is 
approximated by a point dipole. For certain magnetic 
systems their description in terms of a lattice of magnetic 
moments may even be based on the mesoscopic structure 
of the material, especially when a particulate medium is 
described. 



However, our intention is not to describe a particular 
material but to investigate a general model Hamiltonian 
which is 
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where the Sj = fii/fJ> s are three dimensional magnetic 
moments of unit length on a cubic lattice. 

The first sum is the ferromagnetic exchange between 
nearest neighbors with coupling constant J. The second 
sum is the coupling of the spins to an external magnetic 
field -B, the third sum represents a uniaxial anisotropy, 
here, with D e > 0, favoring the z axis as easy axis of the 
system, and the last sum is the dipolar interaction where 
w = ^oMs/(47ra 3 ) describes the strength of the dipole- 
dipole interaction. The are unit vectors pointing from 
lattice site i to j and is the distance between these 
lattice sites in units of the lattice constant a. 

The underlying equation of motion for magnetic mo- 
ments which we consider in the following is the Landau- 
Lifshitz-Gilbert (LLG) equation, 



dt 
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Hi(t)+a(SiXHi(t)) , (2) 



with the gyromagnetic ratio 7 = 1.76 x 10 1:L (Ts) 1 , the 
dimensionless damping constant a (after Gilbert), and 
the internal field H;(f) = -dH/dSi. 

We simulate cylindrical systems being parallel to the 
z-axis with a typical length of 256 lattice sites and dif- 
ferent diameters d. Due to shape as well as crystalline 
anisotropy the equilibrium magnetization is aligned with 
the long axis of the system. However, we start the simula- 
tion with an abrupt, head-to- head domain wall as initial 
configuration, letting the wall relax until a stable state 
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is reached. The distance of the initial wall position from 
the end is approximately 1/3 of the system length. Then 
we switch on the driving magnetic field B along the easy 
axis and wait until a stationary state is reached for some 
time interval in which the velocity v of the wall is con- 
stant while the wall is moving through the central part 
of the wire. We calculate the domain wall velocity from 
the magnetization versus time data, averaged over a pe- 
riod of time where no influence of the finite system size 
on the domain wall can be observed, i. e. until the wall 
approaches the other end of the wire. 

Inspection of the stationary state of the moving do- 
main wall shows that, depending on the ratio uj/J, ei- 
ther transverse domain walls or vortex walls are found. 
Representative spin configurations are shown in Fig. ^ 
The transverse domain wall (left hand side) is observed 
for diameters smaller than the exchange length d ex /a = 
71-^7/(6^(3)) of the system pjj where ((3) « 1.2 is 
Riemann's Zeta-function (see also 1 1 11 for the exchange 
length in continuum theory where 3£(3) is replaced by 
it). Here all spins within cross-sectional planes perpen- 
dicular to the wire axis are parallel so that the system is 
effectively one dimensional. Note, that the spin preces- 
sion leads to a rotation of the spin direction within the 
wall while it is moving. 




FIG. 1: Snapshots of a transverse (left, uj/J = 0.003) and a 
vortex (right, ui/J = 0.2) wall. The diameter d = 8 is kept 
constant while the exchange length of the system is varied. 
Shown is only a part of the system below the current wall 
position. D e / J = 0.05. 

With increasing dipolar interaction, a crossover to a 
vortex wall sets in (right hand side of Fig. ^) which 
is now energetically favorable since the vortex structure 
leads to a flux closure. These findings are in agreement 
with corresponding spin model simulations of therm ally 
activated reversal [T(j and micromagnetic results |12l Il3| 



obtained from simulations of the LLG equation using a 
micromagnetic continuum model. 

In the following we turn to the investigation of the 
influence of the domain wall width and structure on its 
velocity. Fig. |3 compares the dependence of the domain 
wall width A and domain wall velocity v on the strength 
of the dipolar coupling. Here, the domain wall width 
was determined numerically by fitting an tanh-profile to 
the easy axis magnetization of the moving wall in the 
stationary state where the magnetization is averaged over 
cross-sectional planes. However, it should be mentioned 
that for large dipolar interaction in a vortex wall the wall 
profile cannot accurately be described by a simple tanh- 
profile. Note, that even in the limit u> — ► the wall is 
stabilized by the additional crystalline anisotropy D e . 

For a spin chain (d = 1) the domain wall is necessarily 
always planar while for the system with larger diameter 
a crossover to a vortex wall occurs. The crossover can be 
identified as a jump of the domain wall velocity for the 
d = 8 data. Fig. [21 demonstrates that for a transverse 
wall the domain wall velocity is proportional to the wall 
width. For a vortex wall this is at least qualitatively 
the case. Width and velocity of transverse walls decrease 
with increasing dipolar interaction while for vortex walls 
the opposite is true. The crossover itself leads to a jump 
of the wall velocity not the wall width. 
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FIG. 2: Domain wall velocity and domain wall width versus 
dipolar coupling for a spin chain and a cylindrical system, 
respectively. D e /J = 0.05, fi s B/J = 0.05, a = 1. The solid 
lines correspond to Eqs. [I] and [3] error bars are smaller than 
the symbol size. 

For a transverse wall the domain wall velocity is well 
described by an equation derived by Slonczewski as the 
lower limit for one-dimensional domain wall motion ja], 
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Here, A^ is the well-known Bloch wall width 



A R = a, 



J 



2 (D e + 3wC(3)) ' 



(3) 



(4) 
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where, for our case, the denominator D e + 3w£(3) esti- 
mates the effective anisotropy coming from shape as well 
as crystalline contributions (as before in a continuum the- 
ory 3£(3) is replaced by tt). Both equations above are 
drawn in Fig. [5] as solid lines and they agree very well 
with the numerical data for transverse walls. 

In the following we focus on the mobility of vortex 
walls. The crossover from transverse to vortex wall can 
also be observed while varying the diameter of the sys- 
tem keeping the exchange length constant. Since for suf- 
ficiently small driving fields the domain wall velocity is 
proportional to the field in Fig. |3|we directly show the 
domain wall mobility dv/dB versus diameter of the sys- 
tem. Obviously there are two distinct regions with dis- 
tinct wall mobility behavior. For low diameters where the 
transverse wall is found the system behaves effectively 
one-dimensional and in the limit d —>■ 1 the domain wall 
mobility follows Eq. [3] With increasing diameter the 
observed width of the transverse domain wall increases 
little due to the dipolar interaction leading to small de- 
viations from the analytic Slonczewski result assuming 
a Bloch wall width. Nevertheless, we confirmed numeri- 
cally, that Eq. [21 is still valid when the Bloch wall width 
Ab is replaced by the actual (numerically determined) 
width of the transverse wall. 
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stant, namely the Landau-Lifshitz formula 

dv 7 
dB = a ' 



(5) 



where in our case A is the actual domain wall width of 
the vortex wall. The equatioii above is a limit of the more 
general Walker equation 



J 



2{D e + D h sin" 



(6) 



which was derived for sufficiently small driving fields for a 
system with an additional hard-axis anisotropy D^. This 
anisotropy forces the equilibrium magnetization into an 
easy plane. Walkers formula is valid as long as the spin 
motion takes place in one plane which is defined by a 
constant angle <fi to the easy plane of the system, ip is 
given as 



sin cos ( 



y s B 
a2D h ' 



(7) 



where this equation also defines a condition for the va- 
lidity of the walker formula. For a given a there exists a 
maximum field value (or vice versa for a given field a min- 
imum a value) beyond which the spin motion is no longer 
restricted to one plane and instead an irregular preces- 
sional motion starts Note, that the Landau-Lifshitz 
formula is the <fi = limit of the Walker equation, i. c. 
the limit of a strong hard axis anisotropy which forces 
the spin motion into the easy plane. 
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FIG. 3: Domain wall mobility versus diameter of cylindrical 
systems. The model parameters are lu/J — 0.003, D e /J = 
0.1, a = 1. The constant line corresponds to Eq. |31 the 
dashed lines are guides to the eye, error bars are smaller than 
the symbol size. 



Increasing the diameter of the system a crossover from 
transverse to vortex wall is observed with a drastic in- 
crease of the domain wall mobility. As seen in Fig. [3 the 
reason for this effect is not a comparably drastic change 
of the wall width. Instead, as we will discuss in the fol- 
lowing for a vortex wall the domain wall mobility follows 
a law with a different dependence on the damping con- 



FIG. 4: Two different paths for the reversal of a spin. While 
(2) is dominated by precession, as in a transverse wall, the 
path (1) is restricted to one plane as in a vortex wall. 

The equations above were derived for one dimensional 
systems and the question arises why these formula should 
be valid for the motion of a vortex wall with a non- 
uniform spin structure in cross-sectional planes. For a 
qualitative understanding we note that the motion of 
the spins within each spin chain which is parallel to the 
wire axis is indeed restricted to a certain plane passing 
through the spins positions. For a spin chain at the sur- 
face of the cylinder and in the limit of small driving fields 
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these are tangential planes of the cylinder surface. The 
responsible force which keeps the spin motion of each 
chain in this plane is for a vortex wall not a hard axis 
anisotropy — as in the original calculation — but the 
energetical principle which forms the vortex, i. e., the 
combination of exchange and dipolar interaction. Since 
this is the condition under which Walkers formula was 
derived it seems to be plausible that Eq. describes the 
wall mobility in the case of an extended spin system as 
long as the spin motion during the reversal takes place 
in one plane. For a transverse wall, on the other hand, 
the situation is different: the precession of the wall leads 
to the fact that the motion of each single spin consists 
of precession and relaxation with no restriction to one 
single plane. These two different paths for the reversal 
of a spin are sketched in Fig. 

The main difference between Eqs. Eland|Uis the depen- 
dence on damping. This is demonstrated in Fig. [SJwhich 
shows the ratio of the calculated wall mobility and the 
numerically determined domain wall width for two dif- 
ferent strengths of dipolar interaction leading to the two 
different wall shapes. In the high damping limit both 
formulas agree. For a transverse wall the mobility shows 
a maximum at a = 1 and for lower damping the do- 
main wall mobility decreases with decreasing damping 
constant. In the limit a — > only a precession of the do- 
main wall remains without an effective wall motion along 
the wire. 



in our case is played by the combination of exchange and 
dipolar interaction which forms the vortex and forces the 
spin motion into one plane. We would like to stress that 
for experimental systems the low damping limit is more 
relevant. Here, the difference between the two domain 
wall mobilities (reduced to the domain wall width) can 
be extremely large, up to a factor of 1/a 2 . 

For smaller values of a below the critical one the mobil- 
ity decreases again and finally converges to a finite value 
since even for a = the wall can move. In this limit the 
LLG equation conserves the energy of the system, and 
lowering the Zeeman energy leads to an increase of ex- 
change energy leaving an excited spin system behind the 
wall. These observations are also in agreement with the 
calculations of Walker [(| . 
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FIG. 6: Profiles and winding numbers of moving vortex walls 
in a) the high damping limit (a = 1) and b) the low damping 
limit (a = 0). d = 8, u/J = 0.2, D e /J = 0.05, fi s B/J = 0.1 



FIG. 5: Reduced domain wall mobility versus damping con- 
stant. The model parameters are uj/J = 0.01 (transverse 
wall) and ui/ J = 0.7 (vortex wall) respectively, D e /J = 0.05, 
fi s B I J — 0.05, d — 4. Error bars are smaller than the symbol 
size. 



For a vortex wall the domain wall mobility increases 
with decreasing damping constant following a 1/a law 
as long as one is above a critical value a c . As was dis- 
cussed in connection with Eq. [7| this value a c sets the 
limit of pure relaxational spin motion. As was discussed 
before, the role of the hard axis anisotropy Dh in Eq. [7| 



This effect is demonstrated in Fig. Here, profiles 
of the moving walls are shown as well as the so-called 
winding number 
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2nR 



(rotS% dx dy 



which is calculated numerically over all perpendicular 
planes of the wire. The winding number is a measure 
for the existence of vortices, n = 1 means that all spins 
along the boundary are aligned building a ring with per- 
fect flux closure. Fig. Ela) shows the high damping limit 
with a perfect vortex in the center of the wall. In the very 
low-damping limit, Fig. Elb) the situation is much more 
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complicated. Here, behind the moving wall (smaller z) 
an excited spin system is left with vortex-type spin waves 
which are ejected from the moving wall. 

To conclude, in agreement with prior work 
we have found different wall structures for driven domain 
walls in cylindrical systems, transverse and vortex walls, 
depending on the diameter of the system as compared 
to the exchange length of the given material. While for 
vortex walls the domain wall velocity is described by the 
formula from Walker, transverse walls follow a formula 
from Slonczewski. In both cases the domain wall veloc- 
ity is proportional to the domain wall width. The main 
difference is the dependence on the damping constant. 
For small values of the damping constant this difference 
can lead to drastic differences where the velocity of the 
vortex wall is up to a factor of 1/a 2 larger. The reason 



for this difference is probably the fact that each spins mo- 
tion in the case of the vortex wall is completely within 
one single plane as it is the case for the model where 
the Walker formula was derived for, while this is not the 
case for a transverse wall where the precession of the wall 
leads to a three dimensional spin motion. 
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